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Abstract. This chapter is intended to provide a general presentation of 
the atomic and nuclear processes responsible for X-ray line and gamma- 
ray line emission in various astrophysical environments. I consider 
hne production from hot plasmas, from accelerated particle interac- 
^ . tions, from the decay of radioactive nuclei synthesized in stars and 

' from positron annihilation. Spectroscopic properties of these emissions 

, are discussed in the light of the detection capabilities of modern space 

CO ' instruments. 
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1 Introduction 



■ X- and gamma-ray emission lines are valuable signatures of various high-energy 
processes at work in the universe, including heating of astrophysical gas to very 
O ' high temperatures, particle acceleration and nucleosynthesis. Perhaps the most 
^ important realization of pioneering high-energy astronomy was that hot plasmas 
^ are found essentially everywhere in the universe. More and more detailed obser- 
• • vations of thermal X-ray lines provide detailed and often unique information on a 
^ wide variety of astrophysical sites, including stellar environments, supernova ex- 
, plosions, accreting compact objects, interstellar and intergalactic media and active 
\^ ' galactic nuclei. 

Particle acceleration is believed to occur in most of these sites. Our under- 
standing of acceleration mechanisms should greatly benefit from the detection of 
lines produced through atomic and nuclear interactions of energetic particles with 
ambient matter. The sun provides valuable examples of these nonthermal radia- 
tion processes. The study of nuclear de-excitation lines produced in solar flares 
has now become a proper domain of solar physics. 

Gamma-ray lines are also emitted in the decay of radioactive nuclei. Mea- 
surements of gamma-ray activities from cosmic radionuclei testify to the ongoing 
synthesis of chemical elements and their isotopes in the Galaxy and beyond. These 
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observations are now cornerstones for models of nucleosynthesis in novae, super- 
novae and stellar interiors. 

Positron-electron annihilation radiation is relevant to almost all of high-energy 
astrophysics. Positrons can be produced by various pair-creation processes in 
relativistic plasmas, by accelerated particle interactions and by the /S^-decay of 
radioisotopes. Observations of their annihilations offer a unique vision of high- 
energy accreting sources, nucleosynthesis sites and the global structure of the 
Milky Way. 

The aim of this text is to present the physical processes responsible for these 
various line emissions. I also consider emission processes which have not been 
observed with assurance yet: X-ray line production from accelerated particle in- 
teractions and from the decay of cosmic radioactivities, and gamma-ray line emis- 
sion from thermonuclear plasmas. However, given the advent of new X-ray and 
gamma-ray satellites with unprecedented sensitivities and spectral resolutions (see 
Barret and Knodlseder, this volume), I am optimistic that at least some of these 
emissions could be detected in the near future. 

The scope of this review is restricted to emission lines. In particular, I do not 
discuss X-ray absorption lines recently observed from a handful of active galactic 
nuclei, as well as cyclotron absorption features in intense magnetic fields, detected 
from several X-ray binaries. Nor do I treat the high-energy gamma-ray emission 
from pion decay. 

The plan of this review is the following: in section §2,1 present the basic pro- 
cesses of X-ray line emission in both thermally ionized and photoionized plasmas; 
in § 3, I consider nonthermal X-ray line production in interactions of accelerated 
electrons and ions with ambient gas; in § 4, I deal with gamma-ray line emission 
from nuclear collisions and discuss thermonuclear reactions as well as accelerated 
ion interactions; in § 5, I consider both the gamma-ray lines and the X-ray lines 
emitted by the decay of radioactive nuclei synthesized in stars; finally, positron 
annihilation radiation is discussed in § 6. 

2 X-ray line emission from hot plasmas 

Thermal X-ray emission is observed from a great variety of sources in the universe 
(see, e.g.. Astronomy & Astrophysics, volume 365, for an overview of the recent 
XMM-Newton observations). It is produced in various ion-electron plasmas, which 
can be first classified as a function of their optical depth in the X-ray energy range 
(0.1<£;x<100 keV): 

T = t{Ex) = a{Ex) D , (2.1) 

where a{Ex) is the linear absorption coefficient (in c.g.s. units of cm~^) at pho- 
ton energy Ex and D is the typical dimension of the emitting plasma (cm)|^ 
Optically thick sources (t>>1) such as isolated neutron stars essentially emit 



^ a{Ex) is also called the extinction coefficient when it inclu des Compton scattering of the 
X-rays against the free electrons (e.g. Rybicki & Lig htman |l979|). 
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blackbody radiation, with some possible discrete spectral st ructu res due to com- 
plicated radiation transfer effects (e.g. Rajagopal & Romani 1996). Optically thin 
sources include stellar coronas, supernova remnants, superbubbles, the hot inter- 
stellar medium, clusters of galaxies and the intergalactic medium. The emission 
properties of such sources are usually discussed in terms of the coronal model, 
based on the pioneering attempt of Elwert (1952) to explain the X-ray emission of 
the solar corona. In this model, the plasma radiation is essentially due to collisions 
of thermal electrons with ions, but the mechanism for heating the electrons to tem- 
perature Te typically > 10^ K is not specified, i.e. Te is an external parameter of 
the model. 

Different X-ray emitting plasmas are found in astrophysical objects such as 
X-ray binaries and active galactic nuclei, in which a nebular gas of relatively high 
density can be strongly illuminated by the compact X-ray source. The physical 
state and radiation of such photoionized plasmas are generally assumed to be pri- 
marily controled by the emission of the compact source and are usually described 
in terms of the nebular model, which owes its name to the resemblance of these 
objects with planetary nebulae. In this section, I briefly present the X-ray line 
production in the coronal model (§ 2.1) and in the nebular model (§ 2.2) and then 
discuss some diagnostics of plasma parameters based on high-resolution, X-ray line 
spectroscopy (§ 2.3). 



2.1 Thermally ionized plasmas: the coronal model 



The coronal model (see Mewe 
an optically thin, tenuous, thermal plasma. 



199!; for a comprehensive description) applies to 
in which (i) the electrons and the 



ions are relaxed to Maxwellian energy distributions, (ii) the ions are essentially in 
their ground state, i.e. excited state configurations can generally be neglected for 
the calculation of the ionization balance and (iii) each X-ray photon emitted in a 
microscopic process does not further interact with the ambient ions and electrons. 
First order deviations from the coronal model are discussed in Raymond (1988). 



2.1.1 Ionization structure 

The ionization structure of each element of chemical symbol Z and atomic (proton) 
number z is generally obtained by solving a set of z+1 coupled equations: 

# = '^r'C'z+c-) + - r;^(C^+. + , (2.2) 

where 77^ is the ionic fraction of ion Z+' {J2i=oVz = ^)' and az+i (in c.g.s. 

units of cm'^ s~^) are the total ionization (Z+*^Z+(*+^^) and total recombina- 
tion (Z+'^Z+f'^^)) rate coefficients, respectively, and is the electron density 
(cm~^). Note that the above equation is valid only if multiple ionizations can be 
neglected. If the ionization balance is in a steady-state equilibrium [drf^/dt=Q), 
the ionization structure of element Z is simply given by a set of equations con- 
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Fig. 1. Scheme of the dielectronic recombination process: a dielectronic capture (DC) 
into a doubly-excited state a, followed by a stabilizing radiative transition (S) which can 
produce a satellite line (see text). A: Autoionization. 



necting pairs of adjacent ionization states: 



i+i 



(2.3) 



The plasma is said to be at coronal equilibrium when all elements are at ionization 
equilibrium. 

Ionizations are induced by thermal electron collisions. The ionization coeffi- 
cient Cz+i is obtained by averaging the coUisional ionization cross section a^^+i{v) 
over the Maxwell-Boltzmann distribution f{v) of the relative collision velocity, 
which is generally very close to the velocity distribution of the much lighter elec- 
trons: 



C 



with 



fiv) 



-< a 



{v)v >- 



v(Jz+iiv)fiv)dv , 



1/2 



me 



3/2 



exp 



2kn 



(2.4) 



(2.5) 



where rUe and Te are the electron mass and temperature and k is the Boltzmann 
constant. The ionization cross sections are estimated from laboratory measure- 
ments combined wi th ato mic physics theory (in particular distorted wave calcula- 
tions, e.g. Younger 1982 and references therein). 



Apart from direct ionization (Z+'+e ^Z+(*+^'+2e ), ion can be stripped 
of one bound electron through the two-step process of excitation-autoionization 
(Z+'+e"-^(Z+*)*-he"^Z+('+i)-f2e~), which is due to the coUisional excitation 
of an inner-shell electron to a bound level above the first ionization threshold, 
followed by autoionization. The latter process is due to electron-electron interac- 
tions, which can break the compound excited state (Z"*"')* into its two components: 
ion and a free electron. The contribution of excitation-autoionization to 

the total ionization cross section is calculated from the sum over all possible au- 
toionizing levels a of the coUisional excitation cross sections o'z+i i weighted by the 
autoionization probability of each excited state. 
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Fig. 2. Iron ionic fractions at ionization equilibrium as a function of plasma temperature 
(from Arnaud & Raymond 1992| ). 



Recombination of a free electron can proceed either through a radiative free- 
bound transition (Z+^'+^'-l-e"— >Z+*-t-ft,z^) or by a radiationless dielectronic re- 
combination (DR). Charge-transfer reactions are generally negligible for the hot 
plasmas considered in the coronal model, but can be very important for pho- 
toionized plasmas (§ 2.2). Since radiative recombination (RR) is the inverse of 
photoionization, the rates of the two processes are connected through the princi- 
ple of time-reversal invariance. The application of this principle by means of the 
detailed balance theorem leads to the Milne equation (Milne [l924 ) , which is gener- 
ally used to calculate RR rate coefficients from photoionization data (an example 
of the use of the Milne equation is provided in § 2.2.2). 

However, DR is often the dominant recombination process (e.g. Arnaud & 
Rothenflug 1985). It is schematically represented in Figure 1. The first step 
is a dielectronic capture (DC), in which the energy lost by the plasma electron 
is expended in the simultaneous excitation of one of the bound electrons in the 
core of the recombining ion The atomic level a thus produced is said 

to be doubly-excited, because its configuration involves two electrons occupying 
subshells above those which they would occupy in the ground state configuration of 
ion Z"*"'. DC is the inverse of autoionization, so that the DC rates can be calculated 
from the autoionization rates, by applying the principle of detailed balancing. DR 
occurs if, instead of autoionizing, the doubly-excited state decays by a radiative 
transition to a level b below the first ionization limit (Fig. 1). The stabilizing 
decay of the core electron produces a so-called satellite line, as it is slightly shifted, 
usually to the long-wavelength side of the corresponding transition in ion Z~^^^~^^'> , 
due to the electrostatic shielding induced by the spectator outer electron. Satellite 
line emission can provide very powerful diagnostics of plasma parameters (§ 2.3). 

Shown in Figure 2 are the results of Arnaud & Raymond (1992) for the ionic 
fractions of Fe at ionization equilibrium. We see that ions with closed-shell con- 
figurations are more stable than those with partially filled L or M shells. Thus, 
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He-like Fe XXV (Fe+24), Ne-like Fe XVII (Fe+i^) and Ar-like Fe IX (Fe+S), whose 
ground state configurations are Is^, 2s^2p^ and 3s^3p^ respectively, are dominant 
in large temperature ranges, because their ionization and DR rates are relatively 
low compared to those of adjacent ions. 

2.1.2 Line emission 

In the coronal model, the line spectrum is dominated by radiative decays following 
electron impact excitation, plus a smaller contribution of recombination lines. 
Collisional inner-shell ionization can however be important for plasma diagnostics 
based on fine spectroscopic analyses (§ 2.3). Considering only the dominant process 
of collisional excitation, the volume emissivity P^Xi (in units of photons cm^^ s~^) 
of a particular line transition a ^ 6 in ion can be written as 

P^\^ = n.nnazVzS'z^Bab ■ (2.6) 

Here, nn is the H atomic number density (cm^^), is the abundance of element Z 
relative to H, S'§+i (cm^ s^^) is the rate coefficient for electron impact excitation 
of ion Z+* from its ground state to its excited state a and Bat is the radiative 
branching ratio of the transition a b among all possible transitions from level 
a: 

B^, = ^.^^ , (2.7) 

where Aat is the transition probability (s~^) of the spontaneous decay a ^ b (col- 
lisional de-excitation is generally neglected). The excitation rate coefficients 5^+^ 
are calculated by averaging the corresponding cross sections over the Maxwellian 
electron velocity distribution (see eq. 2.4). Semi-empirical fitting formulae have 
been derived for the most significant exciting transitions in all relevant ions (Mewe 



1999| and references therein) . 

The line spectrum is generally dominated by the "allowed" electric dipole (El) 
transitionsQ. However, forbidden lines can be relatively bright if they are pro- 
duced by the decays of excited levels which cannot de-excite by an El transition. 
With the advent of new X-ray instruments with high sensitivity and resolution 
capabilities, such as those aboard the XMM-Newton and Chandra satellites, the 
correct modeling of relatively weak lines has become increasingly important. The 



three main coronal-plasma emission codes, SPEX (Kaastra et al. 1996), CHIANTI 



^ In the Russell-Saunders L5-coupling, which generally holds for the most abundant ions 
(i.e. whos e atn mir, number 2<29), the quantum selection rules of the El transitions are (e.g. 



Sobelman g^79|): AJ=Ja-Jb=0 or ±1, with Ja+Jb>l; AM=Mj^-Mj^=0 or ±1; AL=La-Lb=0 
or ±1, with La+Lf,>l; AS=Sa-Si,=0; na=-nf,; and A£=£-£'=±l. Here, La and Sa are the 
orbital and spin angular momenta of the atomic energy level a, which combine to the total 
angular momentum Ja of magnetic quantum number Mj^ ; na is the parity of level a; and £ and 
£' are the orbital angular momenta of the individual "jumping" electron in the atomic states a 
and b, respectively. Radiative transitions which violate one of these selection rules are said to be 
"forbidden" and are usually weaker. 
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Fig. 3. Thermal X-ray emissions between 0.5 and 10 keV (panel a) and around the Fe 
K line complex (panel b) of an optically thin transient plasma with solar abundances 
and electron temperature Te—SO MK, for three values of the ionization parameter. For 
net=10^'^ cm~'^ s, the plasma has reached the ionization equilibrium. The calculations 
are normalized to the reduced emission measure cm , where D is 

the distance to the source and EM = J neUndV (in cm~^) is the total emission measure 
(which is directly proportional to the X-ray production rate, see eq. 2.6). In each panel, 
the spectra for net=10^^ and 10^"^ cm~^ s are divided by 100 and 10*, respectively. All 
spectra are plotted in energy bins of 10 eV. 



(Landi et al. 1999) and APEC (Smith et al. 2001), now calculate tens of thousands 
of lines in the X-ray energy band. 

As an illustration of the dependence of the X-ray line spectrum on the ioniza- 
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tion structure of a radiating plasma, Figure 3 shows calculated emissions of a tran- 
sient source, that has been instantaneously heated to Te=30xl0^ K at time t—0. 



I used the non-equilibrium ionization (NEI) plasma model (Borkowski et al. 2001 



and references therein), recently implemented in the XSPEC vll spectral software 
package (see Ballet, this volume). Assuming constant electron temperature and 
density, we see from eq. (2.2) that the evolution of the ionization structure can be 
entirely characterized by the product Uet, which is called the ionization parameter. 
The evolving line spectrum is superimposed on a quasi-static continuum, which, 
for Te=30 MK, is mainly produced by the thermal bremsstrahlung of the plasma 
electrons (free-free emission). However, the continuum radiation of plasmas with 
lower temperatures is often dominated by radiative recombinations (Mewe et al. 
198^ ). This emission is considered in § 2.2.2. 



The line emission generally shifts towards higher energies as the emitting ions 
become more and more ionized. For example, we see for n^t—lO^^ cm~^ s in 
panel (a), a prominent line forest at ~1 keV, which mainly consists of Fe L lines, 
i.e. 3^— >2£' electronic transitions in ions Fc XVII-XXIV. This emission is much 
weaker for net=lQ^^ cm~^ s, because in a 30 MK plasma at ionization equilibrium 
the most abundant Fe ion is the He- like Fe XXV (see Fig. 2). The strongest 
line emission is then the Fe Ka line complex (2p— »ls transitions in the Siegbahn 
notation) at '^Q.l keV. Fe K lines are also produced by inner-shell excitation or 
ionization of less ionized Fe, but as shown in panel (6), their energies are slightly 
lower, because the atomic levels get closer as the electrostatic shielding induced 
by the bound electrons increases. For example, the Fe Ka line (which is made of 
the two very close components KQ!i=2p3/2— >lsi/2 a-nd KQ!2=2p]^/2~^lsi/2) is at 
6.40 keV for Fe II-XII, 6.70 keV for Fe XXV and 6.97 keV for Fe XXVI (see panel 
h). Determinations of the precise energy and intensity of various lines thus provide 
direct information on the ionization structure of a coronal plasma. 



2.2 Photoionized plasmas: the nebular model 



The nebular model (e.g. Kallman & McCray 1982 



Liedahl 



199£) is intended to 



describe the physical state and emission of a plasma illuminated by intense X-ray 
radiation. The most important applications of this model are accretion-powered 
X-ray sources: active galactic nuclei. X-ray binaries and cataclysmic variables. In 
these objects, circumsource material is believed to intercept some of the radiation 
emitted from a compact region near the accreting source and to reprocess it into 
line and continuum emission in other portions of the spectrum. The main assump- 
tion of the nebular model is that the temperature and ionization structure of the 
illuminated gas are essentially governed by the interactions of the radiation field 
with the ion-electron plasma. Possible additional heating sources, such as shocks 
or magnetic-field reconnection, are thus generally neglected. 



V. TatischefF: X- and Gamma-Ray Line Emission Processes 9 



2.2.1 Ionization structure and thermal equilibrium 

The temperature and ionization structure of a photoionized plasma are obviously 
intimately connected. For example, recombination of thermal electrons affects 
both the charge state distribution and the plasma cooling rate. Thus, the state of 
the plasma has to be determined by solving a set of coupled equations describing 
the ionization balance together with the thermal equilibrium. Most photoionized- 



plasma models which have been developed so far (see Ferland et al. |1995| and 
references therein) assume that the ionization and energy balances are in a steady- 
state, i.e. that the time scales for variation in the compact source emission and in 
the plasma global parameters are long compared to the characteristic time scales 
for the atomic processes. The ionic fractions of adjacent charge states are then 
connected through a generalization of eq. (2.3), which obviously has to take into 
account the dominant process of photoionization, but also the possible contribution 



of charge-exchange reactions (e.g. Halpern & Grindlay 1980): 



VziCz+^iF^)+neCz+^{T)+K'z+'iT)] - ry^+^ -Kaz+o+D (T)+At+(.+i) (T)] ■ (2.8) 

Here, Cz+'(^£) is the photoionization rate (in units of s~^), which depends on the 
flux of the illuminating radiation (see below) , and i^^+i (T) and K^^i (T) are 
the effective ionization and recombination rates due to charge-transfer reactions 
(also in s^^), which, for a gas of given composition and density, only depend on 
the ionic fractions of H and He, and on the plasma temperature T. Ionization by 
thermal electron impact (of rate coefhcient Cz+i{T), § 2.1.1) is often unimpor- 
tant, because, contrary to coronal plasmas, the temperature of photoionized gas 
may not reach the ionization threshold energy of many ions. Lower temperatures 
can also modify the relative contributions of DR and RR to the total recombina- 
tion rate coefficient a^+it+ijiT) (§ 2.1.1): in a photoionized plasma, radiationless 
dielectronic capture proceeds primarily through core excitations of relatively low 
energies, i.e. essentially transitions of core electrons to the same shell (An=0 ex- 
citations, n being the principal quantum number). Thus, RR generally dominates 
the total recombination rate for H-like and He- like ions, since these species have no 
available An=0 excitations. On the other hand, DR can be t he ma in recombina- 
tion process for ions with partially filled L or M shells (Liedahl 1999| and references 
therein). 

The photoionization rate (^z+^i^e) is obtained by folding the radiation field 
with the corresponding photoionization cross section cr^^i(e): 

CzMF,)= r ^^^^<^lU^)de , (2.9) 



where I^+i is the ionization threshold energy of ion and F^{e,R) is the net 
energy flux at the distance R from the compact source (in c.g.s units of erg s~^ 
cm~^ erg~^). For example, for an isotropic point source of total luminosity Lg and 
spectrum dLg/de=L^s(e), assuming that the source radiation is not significantly 
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attenuated by photoabsorption, the photoioiiization rate is simply given by 

However, photoionization is often enhanced by the emission of one or more 
Auger electrons in addition to the original photoelectron. For example, photoion- 
ization of a K-shell electron in an ion having two or more L-shell electrons leaves 
the resulting ion in an autoionizing state that can decay by ejecting at least one L- 
shell electron. The L-shell vacancy thus produced can further be filled by an other 
radiationless Auger transition if the ion has two or more M-shell electrons. This 
cascade process can lead to the ejection of up to eight Auger electr ons f ollowing 



the creation of a K-shell vacancy in neutral iron (Kaastra & Mewe 1993 ). Thus, 
ionization balance calculations in the framework of the nebular model are gener- 
ally complicated by the possible coupling of several charge states in each element 



(Weisheit 1974). Note that additional ionization can be produced by collisions 



of the fast ejected electrons with neighbouring ions (Halpern & Grindlay 1980 ). 
However, if most of ambient H and He are ionized, these electrons will slow down 
primarily by elastic scattering with thermal electrons, thus depositing their energy 
as heat. 

Charge-transfer reactions can be very important for the ionization structure 
of a photoionized plasma, if a significant abundance of neutral H or He atoms 
(i.e. with neutral fraction typically exceeding '^0.01) coexists with ionized heavy 
elements. They are of the form 

Z+' + A° ^ z+^'-^^ + A+ + AS , (2.11) 

where A=i{ or He and is the energy defect, which can be either a positive or 
a negative quantity depending on the reaction and on the atomic states involved 
in the electron transfer (e.g. Kingdon & Ferland 1999] ). The rate coefficients are 



obtained by integrating the reaction cross sections over the Maxwell-Boltzmann 
distribution of the relative collision velocity. Detailed quantum mechanical calcu- 
lations were performed for few important reactions. For other systems, the cross 
sections can be reasonably estimated by use of the Landau-Zener approximation 



(Kingdon & Ferland 1996 and references therein). 

In the nebular model, the plasma temperature is determined by solving the 
equation of thermal equilibrium: 

UeVe + nH^az^rfz {T^z+^ + ) = 

Z i 

neA^/+nH^az^??z(Af+. +Ag+. +Ag^.) , (2.12) 
z i 

where the left- and the right-hand sides represent the total heating and cooling 
rates (erg cm^'^ s^^) of the thermal bath, respectively. Here, Fg is the heat- 
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ing rate per free electron due to Compton scattering]^ T^+i is the effective pho- 
toionization heating rate, which includes heating by the Auger electrons (ShuU 
1979| ); rg^'.=A'|+°" ^xAE is the heating rate due to exoth ermic {AE>0) charge- 



transfer recombinations or ionizations (Kingdon & Ferland 1999 ); A// is the ther- 
mal bremsstrahlung cooling rate due to free-free emission (see Marcowith, this 
volume); AS+^is the recombination cooling rate due to both RR and DR (Halpern 
& Grindlay 1980 and references therein); A^+i is the cooling rate due to coUisional 



excitation or ionization by thermal electron impacts; and A'^+i is the cooling rate 
due to endothermic {AE<0) charge-transfer reactions. Alternatively, the thermal 
balance can be established by calculating the rate of removal or addition of energy 
to the local radiation field associated with each of the processes affecting atomic 



level populations (Kallman & Bautista 2001) 



2.2.2 Line emission 

It is instructive to compare the line production in photoionized plasmas and in 
thermally ionized plasmas. In the coronal model, the electron population respon- 
sible for the plasma ionization structure is also an efficient source of collisional 
line excitation. On the other hand, a photoionized plasma is generally overionized 
relative to the electron temperature, since the thermal electrons are not the dom- 
inant source of ionization. Collisional excitation is then often unimportant in the 
X-ray energy range and the line emission is preferentially due to recombination 
of relatively low-energy electrons and to photon-induced fluorescence. Therefore, 
whereas in the coronal model, line production in ion is generally proportional 
to its ionic fraction 77^ (eq. 2.6), in the nebular model it is dominated by processes 
involving adjacent charge states: recombination of ion and inner-shell 

photoionization of ion . This difference allows to distinguish between coro- 

nal and nebular gas emission of cosmic origin by observing various line intensity 



ratios from the same element. For example, Liedahl et al. (199C) have shown 
that the {3s^2p) / {3d^2p) line ratios from Fe XVII-XIX can be used to discrim- 
inate between thermally ionized and photoionized plasmas, because the Sd lines 
are preferentially produced by electron impact excitation in a hot coronal plasma 
(5-10 MK, see Fig. 2), whereas the 3s lines are mostly formed by recombination 
in a cooler ('--^0.1 MK) nebular plasma. 

Line emission in the X-ray energy range is dominated by the K lines of C, N, O, 



Compton scattering is generally the dominant heating process of X-ray-emitting photoion- 



ized plasmas. The net rate coef 
cooling can be written as (Ross 



<7X 



at expressing the excess of Compton heating over Compton 



21 



5 nieC^ 



F^{e,R)de 



(2.13) 



where (Tt is the Thomson cross section and c is the speed of light. The first term in the brackets 
represents the heating of relatively cold electrons by scattering with X-rays and includes the first 
order Klein-Nishina correction, and the second term represents the cooling of hot electrons by 
scattering with lower energy photons (see also Marcowith, this volume). 
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Ne, Mg, Si, Ar, Ca, Fe and Ni and the L lines of Fe and Ni. Thus, recombination 
lines are produced in relatively highly ionized ions, having in their ground state at 
least one K-shell vacancy, or one L-shell vacancy for Fe+' and Ni"*"*. On the other 
hand, fluorescent lines being due to radiative transitions from the L or M shell 
following inner-shell photoionization, they are generally associated with colder, 
less ionized plasmas. 

Recombination radiation in overionized plasmas has been observed from sev- 



eral high mass X-ray binaries (e.g. Sako et al. 1999 and references therein). 
It consists of radiative recombination "continua" (see below) and lines produced 
through recombination to excited states followed by radiative cascades. The vol- 
ume emissivity (photons cm~^ s~^) of a given line transition a — > 6 in ion can 
be written as 

P^X, = nenHazriz^^a%+,,+^,Bab , (2.14) 

where a^+t^+i) the effective recombination rate coefficient to the level a (in- 
cluding recombination to higher levels followed by radiative transitions to level a) 
and the other quantities are defined above in eq. (2.6). Radiative recombination 
continuum (RRC) is produced by the first free-bound transition. The monochro- 
matic specific emissivity (photons cm~^ erg~^) from radiative capture of a 
free electron of velocity v is given by 

.ipRRC J 

^±^(e) = nenffaz77^+i«/(z;)af?,.+„ {v)- , (2.15) 

where the energy of the radiated photon is equal to the sum of the ionization 
potential of the recombined ion (Z+*)*, which may or may not be formed in an 
excited state, plus the kinetic energy of the recombining electron, 

e = I{z+')' + ^ruev'^ , (2.16) 

f{v) is the thermal electron velocity distribution (eq. 2.5) and (w) is the 

radiative recombination cross section, which is related to the corresponding pho- 



toionization cross section through the Milne equation (see e.g. Mewe 1999) 



^RR _ . . 



.^+<.+„ (v) = ^iii^ . j-—^ ■ (e) . (2.17) 



Here, (;2'+»=2Jz+>+l is the statistical weight of ion whose total angular 

momentum is J^+i . Combining these equations provides the RRC emissivity as a 
function of the electron temperature: 

dP§^f _ fj ff(z+')» nenHazifz^caf'^+,y{e)e^ ^^^^ I" (e-/(z+.).) 

(2.18) 

The photoionization cross sections are rapidly decreasing functions of photon en- 
ergy from ionization thresholds. For example, they are about proportional to e^"^ 
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for hydrogenic ions (see also Fig. 6a). We see from the above equation that the 
RRC emissivities are steep functions of energy as well. In photoionized plasmas, 
the RRC emissions often appear as relatively narrow features in the spectrum 
(sometimes "line-like" features), with approximate width Ae~fcTe (Liedahl 1999). 
Analyses of RRC shapes can thus allow determination of the electron temperature 
of the emitting plasma. 

Fluorescent line emission has be en observed f rom a number of Galactic X-ray 
binaries (e.g. Ebisa wa et al. 1996| , Sako et al. 199£) and active galactic nuclei 
(Fabian et al. 2000| and references therein). It is produced when an inner-shell 
photoionization is followed by a radiative transition. The volume emissivity of a 
fluorescent line a — > & is thus given by 



b 



(2.19) 



where (P^) is the inner-shell photoionization rate of ion ^\ leading 

to ion in the excited level a' (see eqs. 2.9 and 2.10) and w^+i is the a ^ b 
line fluorescence yield, i.e. the number of photons emitted in the line a b per 
ion in the level a' . Note that the line a ^ & is not necessarily produced in ion 
since the first decay of the level a' can be an Auger transition to an excited 
state a" of ion Auger electron emission can repeat several times before 

the line a — > & is emitted in ion and the calculation of the overall process 

can become cumbersome, since millions of transitions can follow a single inner- 
shell ionization. For example, 5.1 millions of Auger or radiative transitions are 
possible after an Fe I K-shell ionization (Kaastra & Mewe 1993| ). The convention 



is therefore to refer to the inner-shell vacancy produced in ion rather 
than to the excited atomic level a' of ion For example, we speak of the Ka 

fluorescence yield for neutral Fe (for which I will adopt the notation LUpl^=ujp"+a), 
although this quantity strictly describes an atomic property of ion Fe+^. 

The fluorescence yield is a strongly increasing function of ion atomic number 
z. For example, uj^ =0^^"+^^^ varies approximately as ^3.25 f^j. K-shell va- 



cancy produced in neutral atoms up to Fe (Krause 1979). That explains why 
the most prominent fluorescent line observed so far is the Ka line from weakly- 
ionized Fe, since this element has the highest product of Ka fluorescence yield 
(u;^_f^,=0.305±0.001 for Fe I-IX, Kaastra & Mewe | 1993D and cosmic (i.e. so- 
lar) abundance {aFe—3-2xlQ~^, Anders & Grevesse 1989). The line emission at 



~6.4 keV is often accompanied by a sharp absorption edge at 7.1 keV, the K-shell 
electron binding energy in neutral or weakly-ionized Fe: nearly 30% {=ujp^+i) of 
the photons absorbed above this threshold energy are "re-emitted" in the Fe Ka 
line. 

However, if the nebular plasma is optically thick to the ionizing continuum in 
some portions of the spectrum, radiative transfer is generally important and has 
to be taken into account from a geometrical model of the emitting source (e.g. 
Fabian et al. 200C ) . The problem of continuum and line transfer at intermediate 
optical depth is a difficult one and is generally tackled from approximate solutions 
to the full transfer equation (e.g. Kallman & Bautista 2001). It is particularly 
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Fig. 4. Calculated thermal structure (a) and emission (6) of a spherical gas shell of solar 
composition and constant H density (irrespective of its charge state), nif=10^° cm~^, 
illuminated by a steady state, isotropic point source at _R=0. In panel (b), the dashed 
line is the source spectrum, L^f (i?x) = l-45-10''^ xi?x~^ erg s~^ keV~^, and the grey area 
shows the source emission after transmittance through the shell. The total emerging spec- 
trum (solid line) is the sum of the transmitted source emission plus the thermal radiation 
(mostly recombination radiation) and fluorescence emission from the photoionized shell. 



relevant for resonance lines, which correspond to allowed transitions to ground 
ionic levels and can be emitted following radiative recombination or coUisional 
excitation. Indeed, photons emitted in such lines can have a high probability to 
be absorbed by the same transition which emitted them and "re-emitted" at a new 
frequency. This frequency redistribution can repeat many times until the photon 
either escapes the gas or is destroyed by photoabsorption (e.g. Faurobert 1986 ). 

By way of illustrated summary, I show in Figure 4 calculated thermal struc- 
ture and emission of a spherical gas shell photoionized by a central source. I used 
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the XSTAR version 2.1h computer code developed by Kallman et a/.^. We see in 
panel (a) that the temperature is found to decrease from ^10^ K (fcr~8.6 eV) 
in the internal region of the shell to nearly lO"*^ K in the outermost region. Thus, 
collisional excitation is negligible in the X-ray energy range and the line emission 
in panel (&) is mostly due to radiative recombination. The Fe Ka fluorescent 
line emission at ~6.5 keV and the corresponding absorption above 7.1 keV are 
also visible. The irregular continuum emission from the photoionized shell is pro- 
duced by overlapping radiative recombination continua and two-photon decays 
of metastable levels, which cannot decay by a single-photon transition^ (thermal 
bremsstrahlung is negligible in this energy domain). The most prominent RRC in 
panel (&) is the triangular feature just above 54 eV (i.e. in the extreme ultraviolet 
waveband), produced by free electron recombination to the K-shell of fully-ionized 
He. 



2.3 Plasma diagnostics from X-ray line emission 

As already mentioned, simple spectroscopic analyses of the X-ray line emission 
from a hot plasma can provide valuable information, such as the ionization pro- 
cesses (photoionization and/or collisional ionization), the elemental abundances 
and the ionization balance. The recent launch of X-ray instruments with high 
sensitivity and spectral resolution, such as the grating spectrometers aboard the 
XMM-Newton and Chandra satellites, allows more refined diagnostics of plasma 
parameters, among which those based on satellite lines and helium-like line triplets 
may deserve special attention. 

Satellite lines are produced by stabilizing transitions in the process of dielec- 
tronic recombination (§ 2.1.1). The most important of such lines arise from the 
dielectronic capture of a free electron into the L-shell and appear on the long- 
wavelength side of the parent resonance line. For example, the relatively strong 
satellite line j=ls2p'^[^D5/2]^ls'^2p[^P3/2] is at 6.744 A (1.838 keV) for Li-like 
Si XII, whereas the corresponding resonance line w=ls2p[^ Pi]^ls'^[^ Sq] in He- 



like Si XIII is at 6.647 A (1.865 keV; Gabriel [l972| ) . The DR satellite-to-resonance 



line intensity ratios are very sensitive to the electron temperature of a coronal 



" See 

^ The simultaneous emission of two photons is important for H- and 



littp: / /hcasarc. gsfc.nasa.gov/docs/software/xstar/xstar. html 



I He-like ions in the excited 

states n^^^^Lj = 2^ S\/2 ^i^d 2^ Sq, respectively. In this usual notation, n is the principal quan- 
tum number of the outermost occupied shell; L is the total orbital quantum number, which, in 
the Russell-Saunders L5-coupling, is given by the vectorial sum of the individual electron orbital 
angular momenta; S is the total spin angular momentum, obtained by the vectoral coupling of 
the spin angular momentum (s=l/2) of all the electrons (r=2S-|-l is called the multiplicity of the 
level); and J is the total angular momentum, which is the vector sum L+S for LS coupling. The 
orbital quantum number L is usually designated by the symbols S, P, D, F, G. . . (corresponding 
to "sharp" , "principal" , "diffuse" , "fundamental" and the subsequent letters to the alphabetical 
order) for L=0, 1,2,3,4. . .Thus, the level 2'^Si/2 of H-like ions corresponds to the electron being 
in the 2s subshell and the level 2^50 of He-like ions to the two electrons being in the Is and 2s 
subshells with antiparallel spins (5=0). It follows from the quantum selection rules tha t for these 



two levels, the 2s electron cannot decay by a single-photon transition (e.g. Sobelman 1979) 
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plasma. This is fundamentally due to the resonant character of the DR process 
(Fig. 1) compared to the coUisional excitation process: a DR satellite line is pro- 
duced only by free electrons of energy close to the excitation energy of the satellite 
line level (within the small level width), whereas a resonance line can be excited by 
all electrons with any energy above the resonance line energy. This also provides 
a way to evaluate a possible deviation of the electron velocity distribution from 
a Maxwellian distribution. Thus, the existence of a high-energy tail of nonther- 
mal electrons in turbulent plasmas could be identified through an increase of the 
resonance line intensities relative to the DR satellite line intensities (Gabriel & 



Phillips 197c) 



CoUisional inner-shell excitation (IE) can also contribute to the emission of cer- 
tain satellite lines, such as the line (7=ls2s2p[^P3/2]^ls^2s[^ 6*1/2] (at 6.718 A in 
Si XII). Because the simultaneous excitation of two electrons is generally highly 
improbable, some lines (j for example) can only be produced from excited state 
configurations (ls^2p for line j), which are negligible in low-density plasmas. How- 
ever, these excited levels can obtain a significant population at equilibrium for high 
electron density (e.g. at ne~10^^cm^^ for Si XII), thus introducing a density de- 
pendence of satellite lines. 

The intensity ratio of a IE satellite line to its parent resonance line is almost 
independent of temperature, since these two lines are excited by nearly the same 
electron population. However, it does depend on the plasma ionization structure, 
as the intensity of an IE satellite line in ion is proportional to the ionic fraction 
77^1 whereas the intensity of the corresponding resonance line in ion is 
proportional to ?7^^. Thus, satellite lines can also provide diagnostics of the 
ionization state of a transient plasma. 

I now briefly present diagnostics based on the well-known He-like n=2^1 line 
triplet. This system consists of three close lines (Fig. 5): the resonance El line 
w, the forbidden magnetic dipole (Ml) line z and the intercombination line x+y, 
which is made of two unresolvable components that violate the selection rule A5=0 
(see footnote 2) because of magnetic interactions. In Si XIII for example, the lines 
are at 6.684 A , 6.687 A and 6.739 A, for t he x, y and z transitions, respectively. 



As first shown by Gabriel & Jordan ( 1969 ), the emissivity ratio of the forbidden 
to intercombination lines, 

is strongly sensitive to the electron density above a certain value (~5T0^^ cm~'^ 
for Si XIII and ~5-10^ cm""^ for C V). This is due to the coUisional excitation 
of the metastable level 2^S'i to the close levels 2'^Po_i^2, which leads to a decrease 
of the Rz+i ratio with increasing electron density. This density measurement can 
be applied to recombination dominated plasmas as well as to coronal plasmas. 
However, photoexcitation of the 2'^ 5*1 level to the 2^P term may not be negligible 
for photoionized plasmas illuminated by intense ultraviolet radiation (Porquet et 



al. 2001b and references therein). 
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'Pi 




Fig. 5. Scheme of the first energy levels of He-like ions (simplified Gotrian diagram). 
The He- like triplet is formed by the resonance line w, the intercombination line x+y and 
the forbidden line z. Also shown is the continuum-producing, two-photon decay of the 
level 2^5*0. 



Gabriel & Jordan ( 1969 ) also showed that the emissivity ratio 



pz 

Lrz + i — 



(Pi 



p. 



(2.21) 



can be used to diagnose the electron temperature of a coronal plasma, since the 
coUisional excitation rates have not the same temperature dependence for the 
resonance line as for the forbidden and intercombination lines. This ratio can 
also serve as an indication of the validity of the coronal model (e.g. Porquet & 
Dubau 200C). Indeed, a strong resonance line is expected in plasmas dominated 
by coUisional excitation. On the other hand, since radiative recombinations to 
the 2^S and 2^P triplet terms are favored as compared with the 2^Pi level, a 
significant contribution of photoionization leads to a relatively high G^+i ratio 
(see also Liedahl |1999[ ). Porquet & Dubau ( |2000[ ) and Porquet et al. ( pOOlbD have 
recently provided detailed calculations of the line ratios Rz+i and Gz+i for the 
most important He-like ions. 

Satellite lines and He-like triplet lines have proven to be very useful to charac- 
terize laboratory plasmas, as well as coronal plasmas heated by solar flares (e.g. 
Harra-Murnion et al. 1996| ) . The new generation of X-ray spectrometers now al- 
lows to use the same powerful diagnostics for cosmic sources (e.g. Weisskopf et al. 
2002 and references therein; see also Ballet, this volume). 
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3 Nonthermal X-ray line production from accelerated particle interac- 
tions 

Nonthermal electron and ion populations are expected to occur in various astro- 
physical plasmas, when kinetic energy is deposited by acceleration processes into 
the tail of the distributions, at a rate that is sufficiently high to overcome the 
thermal equilibration processes. The impacts of a nonthermal electron tail on the 
ionization balance and X-ray line emission of such "hybrid" plasmas were first 
studied for solar flares and more recently for supernova remnants and clusters of 



galaxies (e.g. Porquet et al. 2001a and references therein). In the following, I 
consider nonthermal X-ray line production in a different context, namely from 
nonthermal particles escaping from the acceleration region and interacting in a 
distinct target medium. It is the most usual situation in laboratory experiments, 
but paradoxically, it has received little attention in astrophysics. 



3.1 Accelerated electron interactions 

Fast electrons interacting in an ambient medium can produce X-ray lines by coUi- 
sional excitation and inner-shell ionization. Nonthermal X-ray line production by 
suprathermal electron interactions has been proposed to account for the observa- 
tion in impulsive solar flares of enhanced Fe Ka line emission at ~6.4 keV that 



is coincident with the hard X-ray burst (Zarro et al. 1992). Recently, Valinia et 



al. (2000) suggested that the interaction of low-energy cosmic-ray electrons with 
interstellar matter could make a significant contribution to the Fe K line complex 
observed in the spectrum of the Galactic ridge X-ray background. 

Here, I assume for simplicity that the interaction region is composed of neutral 
atoms, such that X-ray line emission can result from inner-shell ionization but not 
from collisional excitation. However, the fast electrons can also produce continuum 
hard X-rays via nonthermal bremsstrahlung, which in turn can produce fluores- 
cent line emission via photoionization. Thus, each line can generally result from 
a combination of collisional ionization and photoionization, with relative contri- 
butions depending on the accelerated electron spectrum and on ambient medium 
parameters. I now discuss some fundamental differences between the two atomic 
processes. 

A compilation of cross sections for K-shell ionization is shown in Figure 6a. We 
see that whereas the photoionization cross sections are steep functions of energy 
(they vary approximately as E^^ from ionization thresholds) the cross sections 
for collisional ionization have a much harder energy dependence, in particular, a 



rise in the MeV region caused by relativistic effects (e.g. Hoffmann et al. 1979 ). 
Thus, while fluorescence is essentially produced by photons with energy contained 
in a range of a few keV above the ionization threshold, ultra-relativistic electrons 
can still produce significant X-ray line emission. We also see that the relative 
contributions of collisional ionization and photoionization will generally depend 
on the target element Z. Numerically, the K-shell ionization cross section scales 
approximately as z~^-^ for the electrons, whereas it scales about as z~^-^ for the 
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Fig. 6. (a) Solid curves and data points: cross sections for K-shell ionization of various 
atoms by electron impact. Cross sections calculated from the semi-empirical formula of 



Quarles (1976) are compared with data compiled in Long et al. (1990). Dashed curves: 
cross sections for K-shell photoionization of N and Fe (Verner et al. 1993). (b) Stopping 
range of low-energy electrons in hydrogen. The plotted quantity is Xn/mH, where Xh 
is the stopping range in c.g.s. units of g cm~^ from Berger & Seltzer (1982) and niH is 
the mass of the H-atom. 



photons. 

Propagation makes an other important difference between electron- and photon- 
induced ionization. Figure 6b shows the stopping range of low-energy electrons in 
a neutral H gas. Taking into account ambient He with (n//e/"-//)=0.1 reduces the 
plotted range, measured in H-atoms cm~^, by a factor of ^^1.18. Thus, suprather- 
mal electrons <100 keV injected in a neutral atomic gas of solar composition would 
stop in 4x10^^ H-atoms cm~^. As we will see below, these electrons could produce 
in such a thick target an observable Fe Ka line emission. In comparison, the opti- 
cal depth of >7.1 keV X-rays with respect to photoelectric absorption is Tabs <0. 01 
for iV/f=4xlO^^ cm^^, such that Fc K-shcU ionizing photons would barely interact 
in the target region. 

As a basic introduction to nonthermal processes, let us now calculate the X- 
ray emission produced by the interaction of a steady state, suprathermal electron 
source with a thick target composed of neutral atoms with solar abundances. For 
simplicity, I assume that the nonthermal electron population does not modify the 
state of the ambient medium, i.e. that ionization and heat deposition are negligi- 
ble. I further assume that the typical column density of the interaction region is 
sufficiently low to safely neglect the fluorescent line emission due to photoioniza- 
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Ex (keV) 

Fig. 7. X-ray emission produced by supratliermal electrons witli tlie source spectrum of 
eq. (3.4) interacting in a tliick target composed of a neutral atomic gas of solar abundance. 

The X-ray lines are plotted as delta functions in energy bins of 10 cV. The dashed line 
shows the effect of photoelectric absorption, with a H column density of 4x10^^ cm~'^ 
(see text). 



tion of the ambient atoms by the bremsstrahlung X-rays. The differential X-ray 
production rate (photons s~^ keV"""^) can then be written as 



dQ 
dEx 



(Ex) 
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dNe 
dt 



(Ee) 



riH 



z 



daz 
dEx 



{Ex, El) 



dEi 



idE/dl){E',) 



dEe 



(3.1) 

where (dNe/dt) is the differential injection rate of suprathermal electrons into 
the target region, measured in electrons s~^ keV~^, {daz /dEx) is the differential 
X-ray production cross section for electron interactions with atoms Z (=Z+'^) 
and (dE/dl) is the electron energy loss rate per unit path length in the ambient 
medium, measured in keV cm~^. In this equation, the term in the brackets gives 
the monochromatic X-ray production by an electron of energy Ee as it slows 
down to rest in the ambient medium, which has to be integrated over the energy 
distribution of the nonthermal electrons. 

As H and He are by far the most abundant constituents of the interaction 
region, we have 



dE\ 
dl ) 



ruH 



dE 
dx 



■ aHeiriHe 



H 



dE 
dx 



He. 



(3.2) 



where niH and niHe are the H- and He-atom masses, and {dE/dx)H and {dE/dx)He 
are the electron stopping powers (in units of keV g~^ cm^) in ambient H and He, 
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respectively (Berger & Seltzer 1982| ). Inserting eq. (3.2) into eq. (3.1), we see 



that under the assumption of thick target interactions, the X-ray production rate 
depends on the relative abundances az of the ambient medium constituents, but 
not on the density. 

I consider for the X-ray line production the Ka and K/3 (3p— >ls) lines from 
ambient C, N, O, Ne, Mg, Si, S, Ar, Ca and Fe. The corresponding cross sections 
can be written as 

da^^ 

-^{Ex,E,) = S{Ex - EKi)apE,)u;p , (3.3) 

where Eki is the energy of line Ki (Ka or K/3), 6{Ex — Exi) is Dirac's delta 
function, a^^{Ee) is the cross section for the K-shell ionization of atom Z by an 



electron of energy E^ (Quarles 1976) and ujp is the Ki fluorescence yield for atom 



Z (Firestone 1996t appendix F, and references therein). Note that uj^^=0 for 



z<12 (i.e. Mg), since these atoms do not have 3p electrons in their ground level. 
For the X-ray continuum emission, I take into account electron bremsstrahlung in 
ambient H and He only and calculate the corresponding differential cross sections 



from equation (3BN) in Koch & Motz (1959), with the Elwert correction factor 



which is appropriate at nonrelativistic energies. 

The result is shown in Figure 7 (solid line) for the differential injection rate of 
suprathermal electrons 

— (£;e) = 2.71 • 10*^ X E^^ electrons s"^ keV"\ for 10 < < 100 keV; 
dt 

dN 

—r^(E^ = 0, for E^ < 10 keV or E^ > 100 keV . (3.4) 
dt 

However, the calculated emission is obviously not realistic, since the assumption 
that the target region is sufhcicntly thick to stop <100 keV electrons implies that 
some photoelectric absorption must occur. It generally has to be taken into ac- 
count from a geometrical model of the emitting region. In Figure 7, the dashed line 
spectrum is intended for providing a simple illustration of the effect of photoelec- 
tric absorption. It is obtained by multiplying {dQ/dEx) by exp(— 7V^fcrQf,s(i?x)), 
where (Jabs{Ex) is the absorption cross se ction per H-atom of a medium of cos- 
mic composition (Morrison & McCammon |l983|) and Nh is the mean H column 



density from the site of X-ray production to the edge of the emitting region. I 
adopted for this quantity the stopping range of 100 keV electrons in the ambient 
medium, Nh=4:X10'^^ cm^^. 

The normalization of the source spectrum has been chosen such that the power 
continuously deposited by the electrons in the interaction region is 

E,^[E,)dE, = 1 erg s'^ . (3.5) 

In comparison, the total luminosity of the X-ray emission is 

^.100 kcV iQ 

< / Ex-T^{Ex)dEx = 4.2 x 10"^ erg s"! . (3.6) 

J 0.1 kcV 
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This low radiation yield, Rx=Lx/W, illustrates that X-ray production by supra- 
thermal electrons should generally not be efhcient: in a neutral target region, 
almost all of the electron kinetic energy is preferentially used to ionize ambient 
H and He, whereas if these two species are already ionized, the fast electrons will 
lose most of their energy by collective long-range Coulomb interactions. 

However, the main conclusion of this exercise is that nonthermal bremsstrahlung 
should often be accompanied by observable X-ray lines below 10 keV. In Figure 7, 
the Ka lines from O, Ne, Si and Fe are the four most intense, with calculated 
equivalent widths of 355, 107, 73 and 290 eV, respectively^. A significant line 
emission from elements of relatively low atomic number, such as O and Ne, could 
allow to distinguish X-ray line production by suprathermal electrons from fluores- 
cence emission. This is because in the case of electron impact, the strong decrease 
of the fluorescence yields with decreasing atomic number z is compensated by the 
approximate z~^-^ dependence of the K-shell ionization cross section (Fig. 6a). 



3.2 Accelerated ion interactions 



Low-energy accelerated ions can produce X-rays by a variety of atomic processes 
(Tatischeff et al. 1998 and references therein; Dogiel et al. 1998). A calculated 
nonthermal X-ray spectrum from accelerated ion interactions is shown in Figure 8. 
As for the suprathermal electrons, I have considered a steady state, thick target 
interaction model, in which accelerated ions are injected at a constant rate into 
a neutral ambient medium of solar composition and produce X-rays as they slow 
down to energies below the thresholds of the various reactions. The differential 
X-ray production rate is calculated from an equation similar to eq. (3.1), but 
containing two more summations: one over the elemental abundances of the fast 
ions and the other one over their ionization states. I assumed that the energetic 
particles have the same composition as the Galactic cosmic-ray ions at their ac- 
celeration sources (see Ramaty et al. 1996 , table 1). Their ionic fractions are 



calculated under the assumption that the steady-state nonthermal ion population 
is at t he ion ization equilibrium in the interaction region (see eq. 2.3 and Tatischeff 
et al. 1998 ). I employed the same source spectrum for all ion species. 



dm 

dt 



(3.7) 



which could result from strong shock acceleration in the nonrelativistic energy 
domain (Ramaty et al. 1996). Here, Ei is the kinetic energy of the fast ions, 
measured in MeV nucleon^^, and _Bo is a parameter introduced by Ellison & 



Ramaty (1985) for solar flare acceleration, which is related to the shock size and 
the acceleration time. For both clarity and simplicity, photoelectric absorption in 
the ambient medium is not taken into account. 



^ The equivalent width is the ratio of the Une intensity to the intensity per unit energy interval 
of the underlying continuum at the line energy. 
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Fig. 8. X-ray emission produced by fast ions with cosmic-ray source (CRS) composi- 
tion and the source spectrum of eq. (3.7) with _Eo=30 MeV/nucleon. The calculation 
is normalized to a power of 1 erg s~^ deposited by the accelerated particles in the in- 
teraction region. The total luminosity of the calculated X-ray emission (see eq. 3.6) is 
4.8xl0~^ erg s~^. Note that photoelectric absorption in the ambient medium is not 
taken into account, such that this value is to be considered as an upper limit. IB: inverse 
bremsstrahlung; SEB: secondary electron bremsstrahlung; RRC: radiative recombination 
continuum. 



Continuum X-rays are essentially produced by inverse bremsstrahlung (IB), 
which is the kinematic inverse of normal bremsstrahlung, i.e. the radiation of 
an electron at rest in the moving Coulomb field of a fast ion. Thus, for non- 
relativistic protons of kinetic energy Ep interacting in stationary H, the inverse 
bremsstrahlung cross section is almost identical to the bremsstrahlung cross sec- 
tion for electrons of energy {me/mp)Ep also interacting in stationary H[|. Also 
shown in Figure 8 are the two continuum emissions from the radiative recombina- 
tions of ambient electrons to the K-shell of the fast ions (RRC; Dogiel et al. 1998, 
Tatischeff & Ramaty |1999| ) and fr om the secondary, knock-on electrons that sub- 
sequently produce bremsstrahlung by interacting with the ambient atoms (SEB). 
These two emissions should generally be dominated by the IB production. 

The narrow lines are due to K-shell vacancy production in the ambient atoms 
by the fast ions. In comparison with the similar lines produced by electron impact 
or by photoionization, the lines produced by heavy ion collisions could be shifted 



Baring et al, (2000) have recently shown that in shocked astrophysical plasmas, IB from 
accelerated ions is generally expected to be much lower than bremsstrahlung from suprathermal 
electrons. 
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by several tens of electron-volts, significantly broadened and split up into several 
components, owing to multiple simultaneous ionizations (Garcia et al. 1973)^ 



For example, the neutral Fe Ka line produced by 1.9 MeV/nucleon O impacts is 
blueshifted by ~50 eV in comparison with that produced by proton impacts, and 
has a full width at half-maximum (FWHM) of ~100 eV (see figure 3.55 of Garcia 
et al. L973). Such spectral signatures are within the capabilities of modern X-ray 
instruments. 

The prominent, broad line features in Figure 8 are due to atomic de-excitations 
in the fast ions following electron capture by charge-transfer reactions and coUi- 
sional excitation. We have considered the Ka and K(3 lines from H- and He-like C, 
N, O, Ne, Mg, Si, S and Fe (Tatischeff et al. \L99^ . These lines should generally be 
broad for a steady-state nonthermal population at ionization equilibrium, because 
they are produced by de-excitations of fast moving ion^ For example, the O 
lines in Figure 8 are produced at ~1 MeV/nucleon. The Fe broad lines are emit- 
ted at ~10 MeV/nucleon and merge together to produce a characteristic bump 
in the spectrum, of width =2 keV. However, if the ambient medium is partially 
ionized, since charge exchange could then be considerably reduced, the fast ions 
can significantly slow down before capturing electrons, thus producing narrower 
lines. 

Similar line emission has recently been observed from comets, most probably 
as a result of charge exchange between solar wind particles and cometary neutral 
gazes (e.g. Lisse et al. [2001 ). However, the observed lines are much narrower than 
those of Figure 8, because the solar wind ions remain in high charge states at low 
energies (about 1 keV/nucleon), since they do not reach the ionization equilibrium 
in the interplanetary space. There are as yet no astrophysical observations from 
sources outside of the solar system that would unambiguously indicate an origin 
resulting from relatively low-energy, accelerated ion interactions. Together with 
nonthermal gamma-ray line production (§ 4.2), X-ray line emission from charge 
exchange may probably be one of the most promising signatures of such interac- 
tions and could provide, in particular, one of the best ways of studying Galactic 
cosmic-ray ions at energies below 100 MeV/nucleon. 



4 Gamma-ray line production from nuclear collisions 

Let's now zoom in on the atomic nucleus, from the Angstrom atomic scale to the 
Fermi nuclear scale. Nuclear collisions can lead to gamma-ray line emission in a va- 
riety of ways. These include the direct excitation of nuclear levels, the production 
of excited secondary nuclei, the production of radioactive species decaying into 
excited states of daughter isotopes and the production of neutrons and positrons. 
The latter can originate from the decay of both /^"'"-emitting radionuclei and 7r+ 



* These effects are not taken into account in the present calculations. In Figure 8, the narrow 
lines are simply plotted as delta functions in energy bins of 10 eV. 

^ The present calculations assume isotropic interactions, leading to the maximum possible 
Doppler broadening. 



V. Tatischeff: X- and Gamma-Ray Line Emission Processes 



25 



mesons, as well as from the de-excitation of nuclear levels by e^-e pair emission 



(Kozlovsky et al. 1987). Positron annihilation and the accompanying processes of 
positronium formation and annihilation are fundamentally important in almost all 



of high-energy astrophysics, as first pointed out by Stecker ( |1969| ). It is considered 
in § 6. In the present section, I first discuss the gamma-ray line emission from 
thermonuclear reactions in high-temperature plasmas (§ 4.1), then consider the 
nonthermal gamma-ray line production from accelerated ion interactions (§ 4.2) 
and finally treat the line emission at 2.22 MeV resulting from the radiative capture 
of free neutrons by ambient protons (§ 4.3). 



4.1 Thermonuclear gamma-ray line production 

Thermonuclear reactions are responsible for the nucleosynthesis of the elements in 
stars, but the associated gamma-ray production can not be observed because stars 
are opaque to gamma-rays. Observable thermonuclear gamma-ray line emission 
must be produced in optically thin astrophysical plasmas with ion temperature 
Ti ^ 10^ K {kTi Si 0.086 MeV). Such high-temperature plasmas are believed to 
exist in the vicinity of accreting neutron stars and stellar black holes (Dubus, this 
volume), as well as in central energy sources of active galactic nuclei. In these 
environments, collisions of protons and a-particles (which are the most abundant 
nuclear species) with ^^C and heavier ions can bring nuclei to excited states, which 
can de-excite by a gamma-ray transition. Gamma-ray emission can also be pro- 
duced through radiative nuclear capture, for example by the reaction ^^C(p,7)^^N: 



P 



'■c 



13 



(4.1) 



Recently, Bloemen & Bykov (1997) and Bykov et al. ( 1999| ) proposed that the 
collective gamma-ray line production from an ensemble of unresolved accreting 
neutrons stars could contribute to the diffuse Galactic emission. They were moti- 
vated by COMPTEL observations of the inner Galaxy, which show some evidence 
for an excess emission with a wide latitudinal distribution (^20° FWHM), that 
could cor respo nd to the large-scale height of low-mass X-ray binaries (Bloemen 
& Bykov 1997 ). However, the expected emission from binary sources crucially 
depends on the accretion model. Thus, Bildsten et al. ( 1992 ) predicted very low 
gamma-ray line fluxes from the brightest neutron star X-ray binaries, if the accret- 
ing material is decelerated by Coulomb collisions in the neutron star atmospheres. 
Indeed, these authors showed that the infalling heavy ions would then thermalize 
at higher altitudes in the atmosphere than the accreting protons, such that they 
would be efficiently destroyed by the fast protons, until driftin g to the altitude 
where protons thermalize. On the other hand, Bykov et al. ( 1999 ) considered 
an accretion model in which the infalling ions slow down in a collisionless shock 
produced by the outgoing radiation. In this model, heavy ions can coexist with 
protons and a-particles in the hot postshock plasma for a fraction of second and 
thermonuclear reactions could then produce an observable gamma-ray line emis- 
sion. Bykov et al. suggested that this emission could be detected with INTEGRAL 
from the nearest sub-Eddington sources. 
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Fig. 9. First nuclear levels of (a) ^^C and (b) ^^N (ri/2=9.965 m). Also shown are (a) the 
energy thresholds for the breakup of ^^C into 3 a-particles (7.275 MeV) and into ^^B+p 
(15.957 MeV) and (b) the energy threshold for the breakup of ^^N into ^^C+p (1.944 
MeV). The reaction ^^C(p,7)^^N produces photons of energy (in MeV) £;T.=1.944+-Ecm, 
where Ecm is the total kinetic energy of the ^^C+p interaction in its center-of-mass rest 
frame. 



One of the strongest gamma-ray lines would then be at 4.438 MeV from de- 
excitation of the first excited state of ^^C (Fig. 9a). This level can be populated by 
the inelastic scattering reactions ^^C(p,p')^^C| 439 and ^^C(a,a')^^C4 439, as well 
as by the spallation reactions ^®0(p,pa)^^C| 439 and ^^0(a,2a)^^C| 439!^ Other 
important lines are those resulting from de-excitation of the most abundant heavy 
nuclei, e.g. at 6.129 MeV from de-excitation of ^^O. To illustrate the main char- 
acteristics of thermonuclear gamma-ray line production, I discuss in the following 
the production of the 4.438 MeV line by p+^^C collisions. I also consider the ra- 
diative capture reaction (4.1), which produces a line at 2.365 MeV at the resonance 
energy Ecm=0-'^21 MeV (Fig. 9b). 

I show in Figure 10 the corresponding thermonuclear reaction rates. For the 
^^ C(p ,p'74.438)"^^C reaction, I used the cross section tabulated by Kiener et al. 



(2001) and calculated the thermally- averaged rate from an equation similar to 



eq. (2.4), 



8 \ ^^'^ 1 

<av>^[ TUfTTH / E,m'Jpc{Ec^)e-''-/''^^dE,^ , (4.2) 



The difference of 1 keV between the energy of the gamma-ray line, 4.438 MeV, and the 
energy of the first excited state of ^-^C, 4.439 MeV, is due to the recoil of the ^^C nucleus during 
photon emission. 
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Fig. 10. Maxwellian-averaged reaction rates for p+ C interactions, as a function of the 
ion temperature. Solid curves: productions of the 2.365 and 4.438 MeV hues from the 
reactions ^^C(p,7)^^N and ^^C(p,p'7)^^C, respectively. Dashed curve: destruction of ^^C. 



which is appropriate for an ion plasma of isothermal temperature kTi <^ m^c^ (e.g. 
Clayton 1983). Here, iipc=mpmc / {irip+mc) is the reduced mass of the p+^^C 
system, Ecm is the total kinetic energy of the reaction in its center-of-mass rest 
frame and apC'{Ecm) is the reaction cross section. The rate of the ^^C(p,72.36 5)^'^N 
reaction is obtained from the analytical formula derived by Angulo et al. ( 1999 ) 
for stellar nucleosynthesis purposes, by considering only the contribution of the 
resonance at Ecm=0A21 MeV. We see that the rate of the ^^C(p,p'74.438)"^^C re- 
action several orders of magnitude higher than the one of the radiative capture 
reaction. This is due to the fact that inelastic scattering reactions are mostly 
governed by the strong nuclear force, whereas radiative capture reactions are es- 
sentially electroweak interactions. However, we also see that the 2.365 MeV line 
is produced at lower temperature, because the proton capture reaction occurs at 
lower center-of-mass energy. 

The expected width of the 4.438 MeV line depends on the ve locity distribution 
of the excited ^^C*. Following Higdon & Lingenfelter ( 1977| ), one can neglect 
the p+^^G reaction kinematics and assume that this distribution is essentially the 
same as the Maxwcll-Boltzmann velocity distribution of uncxcited ^^C. The line 
FWHM in keV is then ^9x{Ti/10^ K)i/2, which gives about 90 to 300 keV for 
Ti between 10^° and 10" K. The broadening of the 2.365 MeV line is related to 
the Maxwellian velocity distribution of the center of mass of p+^'^C interactions. 
One has to consider also the relatively large natural width of the first excited state 
of ^^N, r=32 keV (Ajzenberg-Selove 1991). In a first approximation, the 2.365 
MeV line width can then be estimated as FWHM~(r2-t-2.3x lO^'^Ti)!/^, which 
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gives about 35 to 60 keV for T, between 10^ and 10^° K (Fig. 10). In addition, 
lines produced in the vicinity of compact objects can be significantly redshifted, 
but also broadened if the gravitational potential varies over the emission volume 



(Shvartsman 1972) 



Thermonuclear destruction of the heavy species can lead to a significant re- 
duction in the gamma-ray line production (Guessoum & Gould 1988| ; Guessoum 



1989). To illustrate this point, I show in Figure 10 (dashed curve) the reaction 
rate for the breakup of ^^C by proton impact, which I calculated from eq. (4.2) 



using for the destruction cross section the fitting formula of Bildsten et al. ( 1992| , 
eq. C4). This cross section has been estimated by subtracting the 4.438 MeV 
line production cross section from the total p+^'^C inelastic reaction cross section, 
which is justified by the fact that the excited states of -^^C above the 4.439 MeV 
level decay primarily by particle emissio n[3 Let Q be the averaged number of 
4.438 MeV gamma-rays emitted in a thermonuclear plasma by a ^^C nucleus prior 
to destruction. We have (Bildsten et al. |1992D 

Q < , (4.3) 

where <(jv>^ and <av>d are the reaction rates for gamma-ray production and 
^^C destruction, respectively. Thus, we see from Figure 10 that for 7^i3>10^^ 

K, which merely illustrates that destruction, rather than 4.438 MeV line emission, 
is by far the most probable outcome for ^^C nuclei in such very high-temperature 
plasmas. Note that photon- and electron-induced reactions can also destroy heavy 



nuclei. However, Guessoum & Gould (1989) estimated that these processes should 



generally be negligible as compared with breakup by proton collisions. 

As a conclusion, a near future observation of thermonuclear gamma-ray lines 
is uncertain, because the accreting heavy ions could be rapidly destroyed either by 
thermal or nonthermal processes. However, the breakup of ^He and heavier nuclei 
should be accompanied by a copious liberation of neutrons, that may eventually 
recombine with protons to produce a gamma-ray line emission at 2.22 MeV (§ 4.3). 
This line could be produced within the high-temperature plasma near the compact 



object (Bildsten et al. 1993) or in the atmosphere of the companion star (Jean 



& Guessoum 2001). In fact, it may be the best gamma-ray line candidate for 



high-energy accreting sources. 

4.2 Nonthermal gamma-ray line production from accelerated ion interactions 

The interaction of accelerated nuclei (typically >1 MeV/nucleon) with ambient 
matter can produce a wealth of gamma-ray lines with energies ranging from tens 
of keV to about 20 MeV. This nonthermal gamma-ray line emission is often ob- 



served from the Sun during strong solar flares (e.g. Vestrand et al. 199£) and 



it has furnished valuable information on solar ambient abundances, density and 



An exception is tiie 15.11 MeV level, which, because of conservation of isotopic spin, decays 
essentially by 7-ray emission (Fig. 9a). 
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temperature, as well as on accelerated particle composition, spectra and transport 



in the solar atmosphere (e.g. Ramaty & Mandzhavidze 200C; Share & Murphy 



2001). Nonthermal gamma-ray line emission is also expected from various as- 
trophysical sites in which ion acceleration is believed to occur, such as accreting 
compact stellar objects, supernova remnants and active galactic nuclei. Interac- 
tions of Galactic cosmic-ray ions with interstellar matter should produce a diffuse 
emission that might be observed in the near future. These observations would 
then shed new light on the global physical structure of the interstellar medium 
and on the sources of cosmic rays below 100 MeV/nucleon, that are not detected 
near Earth because of the solar modulation. 

The most comprehensive treatment of nuclear de-excitation gamma-ray line 



emission was performed by Ramaty et al. ( 1979 ). The authors took into account 
more than 150 reactions leading to the production of an excited nucleus and eval- 
uated their cross sections from laboratory measurements combined with nuclear 
physics theory. Figure 11 shows a gamma-ray line spectrum which I calculated 
with the corresponding computer codeQ. I have considered the same steady state, 
thick target interaction model as for the nonthermal X-ray emission shown in 
Figure 8, in particular with the same composition (CRS) and source spectrum 
(eq. 3.7) for the accelerated ions. The luminosity of the calculated gamma-ray line 
emission is 



^8 McV jQ 

L^= i E^-^{E^)dE^ = 1.1 X 10"^ erg s~^ , (4.4) 

'o.i McV d-E^ 



which is comparable to the calculated X-ray luminosity (Lx<4.8x 10~^ erg s""'^, see 
Fig. 8). This should generally be the case, except if the accelerated ion spectrum 
is very soft: a high population of low-energy ions (~1 MeV/nucleon) can produce 
nonthermal soft X-rays without producing a significant gamma-ray line emission. 

We see in Figure 11a combination of broad line features from de-excitations of 
fast heavy ions, narrower lines from gamma-ray transitions in heavy nuclei excited 
by energetic protons and a-particles and very narrow lines (dotted lines) from 
spallation-produced long-lived radionuclei, which can come essentially to rest in 
the ambient medium before decaying to an excited state of their daughter nucleus. 
In addition, the prominent broad line feature at ~0.45 MeV is mainly produced 
by the interactions of a-particles with ambient ^He. 

All these lines are superimposed on a continuum- like emission (dashed curve), 
which is due to a large number of unresolved gamma-ray lines, mostly arising 
from cascade transitions in high- lying levels of heavy nuclei. Ramaty et al. ( 1979[ ) 



estimated the intensity an d the energy distribution of this unresolved emission from 



the data of Zobel et al. ( 1968 ). These authors measured the total production of 



gamma-rays of energies >0.7 MeV, in interactions of protons and a-particles with 



The code has been recently updated by taking into account new lal 
improved theoretical e valuations of many cross sections (Kozlovskv et al. 
loaded from the URL: ittp://lheawww. gsfc.nasa.gov/users/ramaty/ViewPubs/ramaty.html 



ory data allowing 
:0n2h. It can be dow n- 
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Fig. 11. Calculated gammarray line emission produced by accelerated ions with the same 
composition (CRS) and source spectrum as in Figure 8, also interacting in a neutral 
ambient medium of solar composition. The calculation is normalized to the same power 
of 1 erg s"'^ deposited by the accelerated particles in the interaction region. The dashed 
curve shows the estimated contribution from all unresolved gamma-ray lines in ions 
heavier than ^^O (see text). The very narrow dotted lines are from the decay of long- 
lived radionuclei. The arrows point to the lines which are mentioned in the text. The 
prompt lines are labelled with the excited nuclei from which the gamma-rays are emitted 
(e.g. ^^C*), whereas the delayed lines are labelled with the parent long-lived radioisotopes 
(e.g. ^®A1). The so-called a— a line feature is mainly produced by a+a fusion reactions 
(see also Fig. 12). 



complex nuclei, at laboratory energies >10 MeV/nucleon. They found that for 

target nuclei heavier than O, the resolvable lines from dc-cxcitation of the low- lying 
nuclear levels can account for only a fraction of the total gamma-ray production. 
In the spectrum of Figure 11, however, the unresolved component produces about 
half of the total gamma-ray emission below 3 MeV. New experiments could allow 
to specify the importance of this emission. 

The broad lines tend to overlap, such that only few well-defined features from 
energetic heavy ions can be distinguished in the gamma-ray spectrum. It is the 
case of the broad line feature from fast ^^C*, centered at 4.4 MeV (Fig. 11) with 
FWHM of ~1.5 MeV, which is produced by interactions of accelerated ^^C and 
with ambient H and He. An other example is furnished by the broad line 
centered at 0.158 MeV with FWHM of ~50 keV, which is due to the de-excitation 
of the first excited state of ^^Co, populated by the reaction in reverse kinematics 
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iH(56Fe,5^Co*)n. 

Intense narrow lines are due to excitations of low-lying nuclear levels in abun- 
dant ambient ions, e.g. at 0.847 MeV from ^^Fe*, 1.37 MeV from ^^Mg*, 1.63 MeV 
from 2°Ne*, 4.44 MeV from ^^C* and 6.13 MeV from ^''O*. Other important nar- 
row lines are those arising from transitions in the spallation products of abundant 
ambient nuclei, e.g. at 0.718 MeV from ^"B*, which is produced by the spallation of 
ambient ^^C and ^^O, and at 0.931 MeV from ^^Fe*, which is essentially produced 
by the reaction ^^Fe{p,pn)^^Fe* . The narrow lines produced in a gaseous ambient 
medium are generally broadened by the recoil velocity of the excited nucleus and 
their FWHM is about 0.5-5% of the transition energy. However, some lines pro- 
duced in interstellar dust grains can be very narrow, because some of the excited 
nuclei can stop in solid materials before emitting gamma-rays^ (Lingenfelter & 
Ramaty 1977). It requires that (i) the mean life time of the excited nuclear level or 
of its radioactive nuclear parent is longer than the slowing down time of the excited 
nucleus in the grain material and (ii) the mean distance from the site of the nu- 
clear interaction to the grain edge is higher than the stopping range of the excited 
nucleus in the grain. The most promising candidates are from the de-excitation 
of the levels of ^^Fe at 0.847 MeV (Ti/2=6.1 ps) and 1.238 MeV (ri/2=640 fs), of 
24Mg at 1.37 MeV (Ti/2=1.35 ps), of ^^Si at 1.78 MeV (Ti/2=475 fs) and of ^^O 
at 6.13 MeV (Ti/2=18.4 ps). Most of the interstellar Fe, Mg and Si, which are 
refractory elements, could be contained in dust grains, whereas about half of the 
interstellar O could be in grains (e.g. Savage & Sembach 1996 ). The detection of 
these very narrow lines, maybe with the INTEGRAL spectrometer (SPI), would 
provide unique information on the composition, size and spatial distribution of the 
interstellar grains. 

Both the broad and the narrow lines could in principle have complex profiles 
that generally depend on the composition, energy and angular distributions of the 
accelerated particles[^, the composition of the ambient medium and the mech- 
anism of the nuclear reactions producing the excited nuclei. In particular, the 
detailed calculation of line shapes requires the knowledge of the angular and en- 
ergy distributions of the recoil nuclei and the angular distribution of the emitted 
gamma-rays. Ramaty et al. (1979) estimated the shapes of the various lines from 
a limited number of available laboratory data together with simplifying theoretical 
assumptions. More recent works have shown that fine spectroscopic analyses of the 
gamma-ray line profiles could give valuable information on the directionality of the 
accelerated particles in the interaction region (Kiener et al. 2001 and references 
therein) . 

In Figure 11, the most intense narrow lines from spallation-produced radionu- 
clei are at 0.478 MeV from the decay of ^Bc (Ti/2=53.3 d), at 0.847 and 1.238 MeV 
from the decay of ^^Co (Ti/2=77.3 d) and at 1.275 MeV from the decay of 22Na 



This is not taken into account in the present calculations, for which I assumed that the 
ambient medium is devoid of dust grains. 

^■^ In the updated nuclear de-excitation line code of Ramaty et al., the accelerated particles 
are assumed to have an isotropic distribution in the interaction region. 
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Fig. 12. Profiles of the gamma-ray line emission at ~0.45 MeV, for 4 values of the 
H density in the interaction region. The total emission (solid curves) is the sum of the 
prompt lines from the reactions *He(a,n7o.429)^Be and ^He(a,p7o.478)^Li (dotted curves) 
and the delayed line at 0.478 MeV from ^Be decay (dashed curves). 



(^1/2=2.6 y). These lines are very narrow only if the density of the ambient 
medium is high enough for the radionuclei to come essentially to rest prior to 
decay. The line at 1.809 MeV from ^^Al decay should generally be very narrow, 
because the time to stop energetic (<100 MeV/nucleon) ^^Al in a medium of 
~1 H-atom cm~'^ is shorter than its half-life (Ti/2=7-4x 10^ y). But ''Be ions, 
which are mostly produced by the reaction *He(a,n)'^Be with recoil energy >2 
MeV/nucleon, stop in a medium with solar abundances in less than 53.3 days only 
if the H density exceeds 10^ cm~'^. 

Figure 12 shows detailed profiles of the gamma-ray emission at ^0.45 MeV for 



lower values of the ambient H density (Tatischeff et al. 2001). The delayed line 
at 0.478 MeV from ''Be decay is superimposed on a prompt emission from the 
reactions ^He(Q;,n7o.429)^Be and ^He(Q;,p7o.478) ''Li- The two prompt lines merge 
(for the assumed isotropic distribution of energetic a-particles) thus producing 
this characteristic emission feature, whose importance in gamma-ray astrophysics 
was first pointed out by Kozlovsky & Ramaty (1974). We see that the width of 
the delayed line decreases as the H density increases, because ^Be ions decay at 
lower energies in a denser propagation region. This provides a potential density 
measurement, which could allow, in particular, to determine if low-energy cosmic- 
ray ions can penetrate into the core of interstellar molecular clouds. 

In summary, interactions of energetic particles with ambient material produce 
a great variety of nuclear de-excitation lines, that could provide a unique tool 
to study the acceleration processes at work in the Galaxy and beyond, and also 
to examine the ultimate nature of the interaction regions. In particular, obser- 
vations of diffuse Galactic emission would permit the first measurements of the 
isotopic composition of the interstellar medium, independent of the ionization and 
molecular states of its constituents. 



V. Tatischeff: X- and Gamma-Ray Line Emission Processes 



33 



B 

10^ 



10^ 



10 



2 3 4 

1 10 10 10 10 

Ej (MeV/nucleon) 

Fig. 13. Neutron production cross sections in units of millibarn=10~'^^ cm~^, as a func- 
tion of tiie proton or a-particle bombarding energy. Here, pp, pa, aa, pHN and aHN 
indicate neutron production in proton-hydrogen, proton-iielium, alpha particle-helium, 
proton-heavy nuclei and alpha particle-heavy nuclei interactions, respectively. The cross 
sections of the three neutron-producing reactions in proton-helium interactions are shown 
separately (dashed curves). The pHN and aHN reaction cross sections take into account 
all isotopes equal to or heavier than ^^C (see text). All the cross sections include the 
neutron multiplicity. 




4.3 The radiative capture of neutrons 

Free neutrons can be produced by the breakup of all isotopes equal to or heav- 
ier than ^H, either by thermonuclear or nonthermal reactions. Their subsequent 
radiative capture by ambient nuclei can produce many specific gamma-ray lines, 
among which the line at 2.22 MeV from the reaction ^H(n,7)^H is expected to be 
the strongest in most astrophysical sites. 

I show in Figure 13 the cross section for important neutron-producing reactions. 



The pHN and aHN reaction cross sections are from Hua & Lingenfeltcr (1987a; 



see also Ramaty et al. 1975). They are defined as 



f^pHN 



Eai<Jp, (4.5) 




O-aHN = ajCTQi (4.6) 
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Fig. 14. Cross sections for the elastic scattering of a neutron on a free ^H-atom and on 
a ^H-atom contained in a water molecule (dashed curve), and for the radiative capture 



reaction '^H(n,7)'^H, as a function of the neutron energy (from Horsley 1966). The radia- 
tive capture cross section satisfies CTcap(-En)=2.44xlO~^//3n (in b) for En<l keV, where 
P„ is the neutron velocity relative to that of light. 



where Ui is the solar abundance of isotope i with respect to ^H, api (resp. aai) is 
the cross section for neutron production in collisions of protons (resp. a-particles) 
with isotope i and the summations are over all isotopes equal to or heavier than 
^^C {J2k afc=l-6x 10~^). The cross sections of the pp, pa and aa reactions are from 



Murphy et al. (1987; see also Meyer 1972). These reactions should generally be 
the dominant neutron-producing processes in interactions of accelerated particles 
of energy >10 MeV/nucleon and in thermonuclear plasmas. However, spallation 
of heavy nuclei plays an important role in solar flares, because the accelerated ions 
have generally a rapidly decreasing energy distribution. 

Shown in Figure 14 are cross sections for n+K interactions. We see that the 
neutrons have a much larger probability to scatter on ambient ^H-atoms than to 
be captured. They lose on average about half of their kinetic energy in each elastic 
scattering. Thus, most of them get thermalized before being possibly captured, if 
they haven't decayed (n^p+e^ or escaped from the interaction region. In a 
medium of H density tih , the mean time for radiative capture is given by 

^ (4.7) 



nHCTcapPnC 

1.4 X 10^^ s 



for En < 1 keV 



(4.8) 



nn (cm -') 

as (Tcap(£'„)=2.44xl0"3"//3„ (in cm^) for E„<1 keV (Fig. 14). This characteristic 
time must be compared with the mean time for neutron decay: 

Tl/2 



Tde 



ln(2) 



(4.9) 



V. Tatischeff: X- and Gamma-Ray Line Emission Processes 



35 



where Ti/2=613.92±0.55 s is the neutron half-life (Groom et al. ^000[ ). Thus, 
significant gamma-ray line emission at 2.22 MeV can only be produced in an 
ambient medium of H density exceeding ~10^^ cm"'^, i.e. in a stellar environment. 
The line should generally be narrow if the only broadening mechanism is thermal. 
If it is the case, its FWHM exceeds the energy resolution of a Germanium detector 
(~2.5 keV at 2.2 MeV) only if the ambient medium temperature is ^ 5x10® K. 

As pointed out by Wang & Ramaty ( 1974 ), ambient ^He can constitute an 
important non-radiative sink for the neutrons. Indeed, the cross section for the 
reaction ^IIe(n,p)^II has the same form at low energies as that of the ^II(n,7)^H 
reaction {a{En) oc l//3n), but is ~1.6xl0'* times larger, which is to be compared 
with the ^H/^He solar abundance ratio, 7.23x10^* (Anders & Grevesse 1989 ). In 
fact, this coincidence has allowed an original determination of the '^He photospheric 
abundance, from observations of the 2.22 MeV line produced in solar flares (Hua 
& Lingenfelter |1987b|). 



The COMPTEL all-sky survey at 2.22 MeV revealed one point-like feature 
(GRO J0332 -87) at the Galacti c coo rdinates (^,6)=(300.5°,-29.6°) (McConnell 
et al. 1997, Schonfelder et al. 2000). The most interesting known source in 



the COMPTEL location error box is RE J0317-853, which is one of the hottest 
(Ri50,000° K) and most highly magnetized («340 MG) white dwarfs, residing at a 
distance of only ~35 pc. Observations with INTEGRAL are expected to provide 
a new insight into the origin of this yet mysterious emission. 



5 Gamma- and X-ray lines from stellar nucleosynthesis 

Observations of gamma-ray lines produced by the decay of radioisotopes synthe- 
sized in stars have long been recognized as a powerful tool for studying the nucle- 



osynthesis processes at work in stellar interiors and explosions (Clayton |1982 and 
references therein; Knodlseder, this volume). However, among the hundreds of 
potential radioactive emitters of gamma-ray lines, only a few of them are suitable 
for astronomical observations. First, nucleosynthesis sites are generally optically 
thick to gamma-rays, such that the appropriate radioisotopes must have a suffi- 
ciently long half- life (typically Tij2>l day) so as to decay only after having left 
the dense production environment. Additionally, their production yields must be 
high enough for the decay activities to be observable by gamma-ray instruments, 
which, in the present state of the gamma-ray line astronomy, practically excludes 
all the radioisotopes beyond the Fe peak. 

Table 1 lists the most important cosmic radioactivities for gamma-ray line ob- 
servations. The astrophysical sites in which these isotopes are synthesized are 
discussed in the revie w paper of Diehl & Timmes (1998). The nuclear data are 
from Firestone ( 1996 , and references therein for the e-/3+ decay branchings), ex- 
cept for the nuclei of mass number A=44, 56 and 57, for which I used the more 
recent data evaluations of Cameron & Singh (1999), Junde (1999) and Bhat ( |1998| ), 
respectively. The half-life of the radioisotopes listed in the Table should be com- 
pared with the characteristic time scale between two nucleosynthesis events which 
produced them in the Galaxy, ~l-2 weeks for novae, ~40 yr for CCSN and ~300 yr 
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Table 1. Star- produced radioisotopes relevant to gamma- ray line astronomy. 



Iso- 
tope 


Production 

sites^ 


Decay 
chain^ 


Half- 
life'^ 


7-ray energy (keV) 
and intensity'* 


^Be 


Nova 


7Be ^ 7Li* 


53.3 d 


478 (0.11) 






56m; ^ ) 56p„* 


U.u f ?J Ll 








56 eiom) 56 


77 9 H 




57Ni 


SNIa, CCSN 


57Ni 57Co* 
57Co ^ 57pg* 


1.48 d 
272 d 


1378 (0.82) 
122 (0.86), 136 (0.11) 




Nova 




2.61 y 


1275 fl) 




SNIa, CCSN 


44Ti 44gg* 


60.0 y 


68 (0.93), 78 (0.96) 






44Sc 44Ca* 


3.97 h 


1157 fl) 


26 Al 


CCSN, WR 
AGE, Nova 


26 Al P'^^^ 26Mg* 


7.4-105 y 


1809 (1) 


60pg 


CCSN 




1.5-10^ y 


59 (0.02f) 






60CO ^ eo^i* 


5.27 y 


1173 (1), 1332 (1) 



Sites which are believed to produce observable 7-ray line emission. Nova: classical 
nova; SNIa: thermonuclear supernova (type la); CCSN: core-collapse supernova; WR: 
Wolf-Rayet star; AGB: asymptotic giant branch star. 

e: orbital electron capture. When an isotope decays by a combination of e and fi^ 
emission, only the most probable decay mode is given, with the corresponding fraction 
in parenthesis. 

Half-lifes of the isotopes decaying by e are for the neutral atoms (see text). 
Number of photons emitted in the 7-ray line per radioactive decay. 
" The underlined 7-ray lines are those for which a positive detection has been reported 
at the time of writing. 

' It is noteworthy that the first excited (isomeric) state of ^^Co at 59 keV decays with 
98% probability by the emission of a conversion electron. 



for SNIa (e.g. Prantzos 1999 and references therein). Thus, we see that we expect 
the decays of the long-lived^ Al and ^°Fe nuclei to produce diffuse gamma-ray line 
emissions, resulting from the superposition of numerous Galactic sources. On the 
other hand, the gamma-ray line activities of ^Be, ^^Ni, ^'^m and '*''Ti should be 
observed in individual transient events. The case of ^^Na is in a way intermediary: 
the line emission at 1.275 MeV could be o bserv ed from nearby individual ONe 
novae (<1.1 kpc with SPI, Hernanz et al. 2001 ), but also from a broad region 
towards the Galactic center, resulting from the cumulated production of few tens 
of active sources. Evidence for this latest unres olved emission has recently been 
reported by the COMPTEL team (lyudin et al. ^002[ ). 

In Table 1, six of the radioisotopes decay mainly or exclusively by orbital 
electron capture (e). As pointed out by Mochizuki et al. (1999) for ''^Ti, one 
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has to be cautious when deducing the production yields of these isotopes from 
the delayed observations of their gamma-ray activities, since their decay rates 
depend on their ionization states. Thus, in the interstellar medium, fully ionized 
e-radioisotopes can only decay by highly improbable nuclear electron capture from 
continuum states, such that they are almost stable as compared with the age of the 
universe^. Mochizuki et al. (1999) suggested that the ionization of ''''Ti behind 
the reverse shock of the Cas A supernova remnant (SNR) could have significantly 
reduced its decay rate, thus leading to a lower mass of ejected ^'^Ti than previously 
inferred from the detected 1.157 MeV line flux. However, a more realistic model 
for the ejecta shows that the effect may not exceed a few percent (Laming 2002). 

The decay of an e-radioisotope can also produce X-ray line emission, when 
the electron vacancy in the daughter atom (primarily in its K-shell) is filled by a 
radiative transition|^. Leising ( pOOl ) has recently examined the exciting possibility 
of detecting cosmic radioactivities in X-rays. Many of them are not gamma-ray 
line emitters, because they decay directly to the ground state of their daughter 
nucleus. The author listed 17 potentially important e-radioisotopes, among which 
the most promising for detection are '^^Fe (Ti/2=2.73 yr), '^''Ti (ri/2=60 yr), ^^Ni 
(Ti/2=7.6xl0'^ yr) and ^^Mn (Ti/2=3.74x 10^ yr). The e decay of ^^Fe could be 
observed through the emission of the Mn Ka line, from very young SNRs such as 
SN1987A. The X-ray line emission which accompanies the ''^Ti decay should be 
hunted for in the two SNRs from which gamma-ray line emission at 1.157 MeV 
has already been observed, Cas A and RX J0852-4622. The detection of this X- 
ray line emission would constitute an interesting complement to the gamma-ray 
line measurements, because it would allow to determine the present ionization 
state of Ti in the two SNRs. Finally, the decays of ^^Ni and ^^Mn, producing 
Ka lines at 6.92 and 5.41 keV, respectively, could be observed from very close 
and reasonably young (i.e. not too extended) SNRs and also as diffuse Galactic 
emissions^. As Leising pointed out, it could be difficult to distinguish the electron 
capture line emission from astrophysical sources which are generally expected to 



^® This is not the case in stellar interiors because of the much higher density of continuum 
electrons. For example, the half-l ife nf ^Be against capture of free electrons in the core of the 
Sun is Ri90 d (Bahcall & Moeller 1961;), which is by coincidence comparable to the measured 
value for neutral atoms. 

The adjustment of the nucleus is generally much faster than the atomic rearrangement, such 
that the atomic transition usually occurs in the daughter element. 

A diffuse emission in the Co Ka line at 6.92 keV is also expected from the decay chain 
13- 60n - 



60 



Fe 



Co {IT: internal transition), since the de-excitation of Co? 



59 kcV 



essentially due to the conversion of a K-shell electron (see Table 1). The expected Galactic flux 
is 

aiK) 



Fx 



Ka 



1.8 X 10 photons s cm 



(5.1) 



1 -I- a 

^ cm~^ is the predicted flux for the 1.173 MeV line from ^"Co decay 
(Timmes & Woosley 1997), a=48 and a(/f)=40 are the total and K-shell conversion coefficients 
for the 59 keV trans ition (King 1993) and a;^^=0.33 is the Ka fluorescence yield for Co. In 
comparison, Leising (2001) predicted that the decay of ^^Ni should produce a 6.92 keV line flux 
from the central steradian of the Galaxy of ~1.7xl0~'' photons s~^ cm~^. 
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be very luminous in X-rays. In particular, detailed models of SNRs are required to 
determine if the expected lines could be unambiguously identified among the other 
thermal and nonthermal X-ray line emissions. However, one can reasonably hope 
that in a few cases, the clear detection of X-ray lines from cosmic radioactivities 
will constitute an important complement to the gamma-ray line studies of stellar 
nucleosynthesis. 



6 Positron annihilation radiation 



Positron-electron annihilation radiation can be produced in a variety of astrophys- 
ical sites. The diffuse Galactic emission at 511 keV is the brightest gamma-ray line 
of cosmic origin. It is produced by the steady-state annihilation of '^1.4-4.3x10'*'^ 
positrons s~* (depending on the spatial distribution of the emission, Kinzer et 
al. 2001). A large fraction of these positrons could result from the decay of 
^^Co and '^^Sc (synthesized as ^^Ni and '^^Ti in supernovae, see Table 1), with a 
smaller contrib ution from the /?+ d ecay of ^®A1 (Chan & Lingenfelter 1993; Der- 
mer & Murphy 2001). Milne et al. (2002) have recently estimated that these three 
isotopes could supply 30-50% of the Galactic positrons. The decay of 7r+ and ra- 
dioisotopes produced by cosmic-ray interactions in the interstellar medium could 
make an additional contribution of a few percent (e.g. Ramaty & Lingenfelter 



19791) . 



Positrons are also produced in classical novae, from the decay of radioisotopes 
synthesized during the thermonuclear runaway. The main /3"''-radionuclei synthe- 
sized in novae are (ri/2=9.965 min), ^^F (ri/2=109.8 min), 22Na (ri/2=2.61 h) 
and 2^A1 (Ti/2=7.4x lO'^ y). The decay of the latter two isotopes could inject a 
small but non-negligible amount of positrons in the interstellar medium, but ^'^N 
and ^^F have too short half-lifes, so that most of the positrons should annihilate in 
the relatively dense material ejected by the nova explosion. The resulting gamma- 
ray emission could be detected with SPI from relatively nearby novae during the 
first day after the explosior 



2001| ). 



(Hcrnanz et al. 

Positron production can also occur in accreting compact objects, from the pair 
production reaction 7-1-7 — > +e+, from the materialization of photons in strong 
magnetic fields 7-1-5 e~+e+ + il and from other less important processes 
such as 7 -|- AT — > + e+ + A, where A = or nucleus. A significant fraction 
of these positrons could annihilate near the compact objects, thus producing a 
localized and possibly time-variable emission. The annihilation radiation from 
relativistic plasmas is discussed by Marcowith (this volume). Here, I consider the 
emission produced by the injection of nonthermal positrons in various astrophysical 
environments. 

Important processes leading to gamma-ray emission from positron annihilation 
are shown in Figure 15 (see Guessoum et al. 1997 for discussion of the various 



The detectability of this emission crucially depends on the uncertain rates of the reactions 
^^F(p,7)l^Nc_and ^^F{p,a)^^0, which are responsible for the thermonuclear destruction of 
(Coc et al. 



200' 



V. Tatischeff: X- and Gamma-Ray Line Emission Processes 



39 










Positroniiim 


Ortho- 


Para- 


1 

3Y 


1 

27 




Fig. 15. Main processes leading to gamma-ray emission from positron annihilation 



(adapted from Guessoum et al. 1991, fig. 1) 



processes and update of the cross sections) . Positrons produced by the decay of 
radioactive nuclei have initial kinetic energies typically ranging from several hun- 
dred keV to few MeV. They slow down mainly by Coulomb interactions with free 
electrons and by excitations and ionizations of neutral atoms, and their character- 
istic slowing-down time is ^10^ /n years, where n is the ambient medium density 
in units of cm~^ (Guessoum et al. 1991 ). Positrons produced by other processes 
generally have higher initial energies. Having lost the bulk of its kinetic energy, 
a positron can either form a positronium (Ps) atom in flight by charge exc hange 
with an atom or a molecule (at ^ 50 eV in neutral H or H2, Bussard et al. 1979 ) 
or become thcrmalized with the free electrons. Thermal positrons either annihilate 
directly with free or bound electrons or form Ps atoms through radiative recom- 
bination and charge exchange (Fig. 15). As first pointed out by Zurek ( |1985| ), 
interstellar dust could play a significant role in annihilation of thermal positrons 
despite its low abundance, since its projected surface area (~10^^^ cm^ per Galac- 
tic H-atom) can be larger than the cross sections for the other processes leading to 
positron annihilation. Collisions of thermal positrons with dust grains can lead to 
the reflection of the positrons from the grains, to the annihilation of the positrons 
within the grains and to the formation and subsequent escape of Ps atoms back 
to the gas. The contribution of these processes to gamma-ray emission crucially 
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depends on the grain properties, including their compositions, sizes, charges and 
abundances in the ambient medium. 

The positronium atom, which is the simplest hydrogen-like system, is prefer- 
entially formed in the state of principal quantum number n=l, either as a triplet 
i^Si) in which the spins of the electron and positron are parallel, or as a sin- 
glet (^Sq) in which they are antiparalleQ These states are called, respectively, 
orthopositronium and parapositronium, by analogy with the spectral designa- 
tions for hydrogen. Parapositronium annihilates with a lifetime Tp=1.25xl0^^*' s 
into two photons of 511 keV in its rest frame, but orthopositronium decays with 
ro=1.4xl0^^ s into three photons, which form a characteristic continuum at en- 



ergies below 511 keV (Ore & Powell 1949 ). Ps atoms can be quenched before 

annihilation if the ambient medium density is ^ 10"'^^ cm'^. 

Let If,+ be the total number of positrons which annihilate in a given source, and 
and 1^^ the number of photons emitted from this source in the 511 keV line 

and in the ortho-Ps continuum, respectively. Since the probability for formation of 

the lepton atom in the triplet state is three times the probability for its formation 

in the singlet state, we have 

- Ie+ X [2 X (1/4) X fp, + 2 X (1 - /p,)] (6.1) 

and = /e+ X 3 X (3/4) x fp, , (6.2) 

where fps is the fraction of positrons that annihilate via Ps formation. Combining 
these two equations gives fps in terms of the observable quantity [12-^ / I^-y): 

2 

^ 2.25(/2,//3,) + 1.5 ■ ^^-^^ 

To illustrate how an observation of the 511 keV line profile and the determina- 
tion of fps can provide unique information on the annihilation environment, I show 
in Figure 16 calculated spectra of the radiation produced by positron annihilation 
in the four phases of the interstellar medium. I used equation (9) of Guessoum 



et al. (1991) with the rates of the various processes and the corresponding line 
widths given in that paper, except for the widths of the lines produced by radiative 
recombination and direct annihilation of thermal positrons, for which I used the 



more accurate results of Crannell et al. (1976; see Guessoum et al. 1997). For 
simplicity, I neglected the processes involving interstellar dust. 

We see in panel (a) that in the cold phase, the line is relatively broad and 
the fraction of positrons which annihilate via Ps formation is high. Indeed, in 
this neutral environment, most of the nonthermal positrons form Ps in flight by 
charge exchange and the immediate annihilation of para-Ps produces a broad line 



(6.4 keV FWHM, Brown et al. 1984). On this component is superimposed a 



The formation of a Ps atom in an excited state can produce a Lya line at 2430 A and 
other hnes in the UV, IR and radio bands, that could prov ide a n other way to observe positron 
annihilation from astrophysical sources (e.g. Wallyn et al. L996). 
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Fig. 16. Calculated gamma-ray spectra of the positron annihilation radiation produced 
in the various phases of the interstellar medium. Processes involving interstellar dust are 
not taken into account. The calculations are normalized to one positron annihilation in 
each phase, fps : fraction of positrons that annihilate through Ps formation. 



narrower line due to the direct annihilation of 10% of the positrons with bound 
electrons. In the warm neutral phase (of temperature T=8000 K and ioniza- 
tion fraction C="-e- /("-e- )=0-15i McKee & Ostriker 1977), the fraction of 
positrons forming Ps in flight is reduced to 22%. Most of the positrons get ther- 
malizcd and form Ps by charge exchange just above the energy threshold (6.8 eV 
in the center of mass fo r e++ H^Ps+p). The annihilation line has a FWHM of 
1.5 keV (Bussard et al. |l979|) . In the warm ionized phase (T=8000 K, C=0.68), 
the 511 keV line is even narrower (1.1 keV FWHM), because a significant fraction 
of thermal positrons form Ps through radiative recombination, which does not 
have an energy threshold. Finally, in the hot phase of the interstellar medium 
(T=4.5xl0^ K), all the positrons get thermalized and the annihilation line should 
be broad («7.4 keV FWHM) due to the thermal moti on of the the center of mass 
of the e+-e~ pairs which annihilate. Crannell et al. ( 1976| ) estimated the width 
of the 511 keV line from both para-Ps annihilation following radiative recombina- 
tion and direct annihilation with free electrons as FWHM=llx (T/IO^ K)^/^ keV. 
However, the line width could be significantly reduced if annihilation on interstel- 
lar dust is important (Guessoum et al. 1991). The results shown in Figure 16 
are in good agreement with the recent calculations of Dermer & Murphy ( ^001 ) , 
except for the cold phase, for which they found the narrow line component from 
positron annihilation with bound electrons to be negligible. 



42 



Title : will be set by the publisher 



Harris et al. (1998) have measured, with the TGRS instrument (Germanium 
detector) on the WIND spacecraft, the width of the 511 keV line from the Galactic 
center region to be narrow, FWHM=1.81±0.54±0.14 keV, and the Ps fraction to 
be large, /ps=0.94±0.04, in good agreement with the most recent OSSE result: 
/ps=0.93±0.04 (Kinzer et al. 2001). The narrow line width suggests that most of 
the positrons can not penetrate deeply into molecular clouds and that annihilation 
in the hot interstellar phase is not predominant, unless dust grains sprinkled into 
the hot plasma are important sites of annihilation. However, the Ps fraction would 
then be very low (<0.5, Guessoum et al. 1991), in conflict with the measured val- 
ues. Thus, these observations favor a scenario in which the Galactic positrons 
annihilate predominantly in the warm, partially ionized interstellar medium. The 
spectroscopic-imaging capability of INTEGRAL should allow a better determi- 
nation of the nature and the spatial distribution of the emitting regions. This 
would shed new light on the sources of the positrons and on the thermodynamic 
properties of the annihilation media. 



I would like to thank Jean Ballet, Jiirgen Kiener, David Lunney and Alexandre Marcowith 
for critical reading of the manuscript. 
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